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ABSTRACT 

Finite cut-off effects strongly influence the thermodynamics of lattice regularized 
QCD at high temperature in the standard Wilson formulation. We analyze the 
reduction of finite cut-off effects in formulations of the thermodynamics of SU(N) 
gauge theories with three different 0(a 2 ) and 0(a 4 ) improved actions. We calculate 
the energy density and pressure on finite lattices in leading order weak coupling 
perturbation theory (T — > oo) and perform Monte Carlo simulations with improved 
SU(3) actions at non-zero g 2 . Already on lattices with temporal extent N T = 4 
we find a strong reduction of finite cut-off effects in the high temperature limit, 
which persists also down to temperatures a few times the deconfinement transition 
temperature. 
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1 Introduction 



With increasing accuracy of lattice calculations it also becomes increasingly impor- 
tant to control and, as far as possible, eliminate the systematic errors introduced 
through the finite lattice cut-off in numerical calculations. Various procedures to 
achieve this goal without increasing the necessary computational effort drastically 
have been suggested in the past, including the early Symanzik improvement pro- 
gramme 0, renormalization group improved actions @ as well as the recently con- 
structed perfect actions 0, [|]. The resulting actions are not designed to lead to 
a better infrared behaviour of lattice regularized field theories. It thus may not 
have come as a surprise, that the application of these actions in lattice calcula- 
tions, which aimed at a determination of long-distance properties of asymptotically 
free field theories, did show little advantages over the standard one-plaquette ac- 
tion originally proposed by Wilson [Q. In particular, this is the case for studies 
of the finite temperature deconfinement transition in SU (N) gauge theories with 
Symanzik-improved actions [^|, |7|. On the other hand the improved actions are ex- 
pected to lead to a better ultraviolet behaviour of the theory. This has, for instance, 
been observed in studies of topological properties f| or the short distance part of 
the heavy quark potential [f|]. In these cases the unwanted short distance lattice 
artifacts could successfully be eliminated. 

The high temperature behaviour of QCD is close to that of an ideal gas. Bulk 
thermodynamic quantities are therefore dominated by contributions from large mo- 
menta. These, however, are most strongly influenced by finite cut-off effects. One 
thus may expect that improved actions will be particularly useful for the analysis 
of bulk thermodynamic quantities. In fact, it is well known that in the standard 
Wilson formulation of lattice QCD the cut-off effects lead to strong modifications of 
the high temperature limit of quantities like, e.g. the energy density and pressure 
H [l(J. These cut-off effects are 0((aT) 2 ) in the pure gauge sector. Calculations 



on lattices with temporal extent N T = 1/aT have recently been performed for the 
SU(2) 0] and SU(3) gauge theories. They clearly show the influence of a 
finite ultraviolet cut-off on the behaviour of bulk thermodynamic quantities and 
their extrapolation to the continuum limit. In the standard Wilson formulation || 
lattices with temporal extent N T >8 are needed in order to reduce deviations of, 
eg. e/T 4 , from the continuum extrapolation below a few percent. The interesting 
bulk thermodynamic quantities like energy density and pressure all have dimension 
four. In lattice calculations they are determined from operators, whose expectation 
values are proportional to N~ A = (aT) 4 . It thus rapidly becomes difficult to cal- 
culate these operators with reasonable statistical significance. In fact, the effort to 
calculate them with a given accuracy on lattices of size N% x N T , increases at least 
like iV 11 , if one also keeps the physical size of the thermodynamic system constant 
(N a /N T = const.). It therefore should be clear that a huge improvement is already 
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achieved, if one can perform calculations on lattices with, say N T = 4, with cut-off 
distortions similar in magnitude to what one reaches in calculations with the Wilson 
action on lattices with twice that temporal extent. 

It is the purpose of this paper to quantify the ultraviolet cut-off effects that can 
be expected to be present in calculations of thermodynamic quantities with improved 
actions and to examine the relevance of improved actions for the calculation of the 
equation of state in QCD at high temperature. We will discuss various improved 
actions and calculate finite cut-off corrections to the high temperature limit of the 
energy density in leading order perturbation theory. We also will present some 
numerical results for the pressure of a SU (3) gluon gas calculated with an 0((aT) 2 ) 
improved action. We are going to discuss tree level improved actions (g 2 = 0). As 
will become clear in the following this leads already to a significant improvement of 
the high temperature behaviour of bulk thermodynamic quantities in SU(N) gauge 
theories. There is, however, no fundamental problem which would prohibit the 
extension of the present considerations to 0(g 2 ) improved actions and we will present 
evidence that this is, in fact, needed for temperatures close to the deconfinement 
transition temperature. 

We will present in the next section two extensions of the standard one-plaquette 
Wilson action leading to an 0((aT) 2 ) improvement of thermodynamic observables 
as well as a specific choice of an action that leads to an improvement at 0((aT) 4 ). In 
section 3 we will discuss the perturbative high temperature limit of these improved 
actions. A first exploratory numerical analysis of these actions and a comparison 
with the standard one-plaquette Wilson action is presented in section 4. In section 
5 we give our conclusions. Various details of our perturbative calculations are given 
in an Appendix. 



In order to eliminate the 0(a 2 ) and higher order corrections to the lattice version 
of the Euclidean action one can add appropriately chosen larger loops to the basic 
4-link plaquette term appearing in the standard one-plaquette Wilson action. We 
will, in particular, discuss here simple extensions of the one-plaquette action, which 
only involve larger planar loops, i.e. we will consider the generalized Wilson actions, 



2 Improved SU(N) Actions 



oo 



s 1 = 



E EE *knw%(z) 



x,v>n 1=1 k=l 



E s^{x) 



(2-1) 
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where W^' l u denotes a symmetrized combination of k x I Wilson loops in the (/i, v)- 
plane of the lattice, 

Wfttx) = 1 - ^ (Re TrL<»L® k ^^LW+ + (* <-> /)) . (2.2) 

Here we have introduced the short hand notation for long links, M = Ilf=o U x+ j^ 
and x = (nx, n 2 , n 3 , n 4 ) denotes the sites on an asymmetric lattice of size N% x iV T . 

With a suitable choice of the coefficients a^i it can be achieved that the gener- 
alized Wilson actions reproduce the continuum Euclidean Yang-Mills Lagrangian, 
C = — \F^ V F^ V , up to some order 0(a 2n ) |[J. The standard one-plaquette Wilson ac- 
tion, S^ 1 ' 1 ), with ai 5 i = 1 and a^j = for all (k, I) ^ (1, 1) receives 0(a 2 ) corrections 
in the naive continuum limit. Expanding the link variables, U x ^ = expiigaA^x)), 
in powers of a one finds 

= —^(f^F^ + ^F^{dl + d 2 u )F^ 

+0(a 4 )) . (2.3) 

By adding an additional Wilson loop to the action one can achieve that corrections 
start only at 0(a A ). In particular we will consider here actions obtained by adding 
a planar 6-link and 8-link loop, respectively. The non-vanishing coefficients in these 
cases are, 



J = (1,2): <zi,i = - , a 1)2 = -~ 

6 O 

I =(2, 2): a ^ = \ > a2 ' 2 = ~^8 (2 ' 4) 

The action S^ 1,2 ) is a specific choice of the 6-link improved actions originally proposed 
by Symanzik The action S*( 2 ' 2 ) has recently also been discussed in the context 
of NRQCD calculations [0, jlB] . Its generalization to larger quadratic loops gives a 
straightforward procedure to eliminate also higher order cut-off effects. Already in 
0(a 4 ) there exist two independent operators of dimension eight, which contribute 
to the finite cut-off effects. In general, one thus needs three independent loops of 
length six and eight in order to eliminate all cut-off effects proportional to a 2 and 
a 4 , respectively. In an action constructed only from quadratic loops the different 
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operators contribute, however, with the same relative weight. In that case it thus is 
sufficient to add only two additional, quadratic loops to obtain an 0(a 4 ) improved 
action. We will consider here the case where we add 2x2 and 3x3 Wilson loops 
to the plaquette termF|. In this case corrections to the continuum limit start only at 
0(a 6 ). The non- vanishing coefficients for the action S^ 3 ' 3 ) are given by, 

3 3 1 

I = (3,3) : ai,i = - , a 2 , 2 = — , a 3i3 = — (2.5) 



In the following we will discuss the behaviour of thermodynamic quantities in 
the high temperature limit, which are obtained from the partition function, 

Z= [U dU ^ e ~ PSl > ( 2 - 6 ) 

defined on lattices of size A*" 3 x N T . Here (3 = 2N/g 2 is the bare gauge coupling. 
Thermodynamic quantities like the energy density or the pressure are obtained as 
derivatives of the logarithm of the partition function with respect to the temper- 
ature, T = l/N T a, and the volume, V = (N a a) 3 . They are thus represented by 
expectation values of certain parts of the action, i.e. the operators used to evaluate 
thermodynamic quantities are improved to the same order as the action. In the 
following we will show explicitly that the 0(a n ) improvement of the action will lead 
to an 0((aT) n ) improvement of thermodynamic observables. 



3 Perturbative High Temperature Limit 
3.1 The generalized SU(N) Wilson Actions 

A suitable quantity for the study of the influence of finite cut-off effects on the 
thermodynamics of SU(N) gauge theories is the leading high temperature behaviour 
of the energy density or equivalently the pressure. In the limit T — > oo the energy 
density approaches that of a massless, non- interacting gluon gas, e/T 4 — > €sb/T 4 = 
(N 2 — l)7r 2 /15. Corrections to this are of 0(g 2 ). Numerical investigations of the 
high temperature phase of SU(N) gauge theories show that the energy density 
rapidly approaches this high temperature limit. More precisely, on lattices with finite 
temporal extent N T , it approaches a limiting value, which differs from €sb/T a due 

a This procedure can easily be extended to higher orders. The resulting actions are, however, 
non-local in the sense that the magnitude of the coefficients of a k x k loop asymptotically only 
decreases powerlike, \au,k \ ~ k~ A . 
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to ultraviolet cut-off effects. At temperatures a few times the critical temperature 
deviations from this limiting value are generally of the same order of magnitude 
as the cut-off distortion effects. In order to quantify the approach to the high 
temperature ideal gas limit it thus is important to reach a good approximation of 
this limit already on finite lattices. The cut-off effects due to finite values of N T 
can be studied in leading order perturbation theory. In fact, in the case of the 
one-plaquette Wilson action this is known for quite some time || |10 . 



For the analysis of the high temperature ideal gas limit it is sufficient to evaluate 
perturbatively the difference between expectation values of space- and timelike parts 
of the action, 



Ti ^h = ^K{{sl)-{s- T )) + o(i-) 



\(!P-l)N$(Sl^-Xi/-')i Our , . 13.1; 



where S^ 2 ^ and S^ 2 ^ denote the 0(g 2 ) expansion coefficients of the expectation 
values of the space- and timelike parts of the action, 



1 N 2 — 1 

<#> - ^( E 4 4 (,)>^ 2 y^) + 0(^) . (3.2) 



All other terms appearing in the definition of the energy density are 0(g 2 ) and thus 
do not contribute in leading order perturbation theory JET 



In the following we will discuss the perturbative calculation of the expansion 
coefficients S^ 2 ' and S^ 2 ' on finite lattices. Thereby we will also review some of the 
results obtained in perturbative calculations with the Wilson one-plaquette action 
on finite, asymmetric lattices || Our notation closely follows ||10|| . Perturbative 
results for the Symanzik action, S^' 2 \ on infinite lattices have also been obtained 



in 11 . 



The improved Wilson actions are defined in terms of the link variables^, U XtlM = 
exp (igA^x)). They can be expanded in powers of g 2 and represented in momentum 



3 In the following the cut-off a has been set to unity. 
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space. Sums in Fourier space we denote by 



| = , (3-3) 



where p M = (27r/iV -)n /i , ra M = 0, 1,..., iV CT — 1 for \i = 1, 2, 3 and = (27T / 'N T )ri4, 
7i4 = 0, 1,..., N T — 1. The lattice gauge fields are defined on the middle of a link, 

MP) = E e_W " ip " /2 M* + M*) > ( 3 - 4 ) 



with A M = A£A a , a = 1, (N 2 - 1), and the normalization 2TrA a A fe = 6 ab . We also 
will use the short hand notation = sin(p M /2). 

For the evaluation of the leading order term in the energy density on finite 
lattices one only needs to keep the quadratic part in the perturbative expansion of 
the action, (3S 1 , and introduce a gauge fixing term. We combine both in the form 

So = -~ / E^C-p^CpKCp) > ( 3 - 5 ) 



with the inverse propagator 

A M) „(p) = G^(p) + ig^gM ■ (3.6) 

The first term in Eq. (|3.6|) arises from the expansion of the action and the second 
term defines the covariant gauge fixing. It is convenient to separate in G^ v a diagonal 
part, 

G„J P ) = DJp)5„„ -E.jp) , (3.7) 



with 



oo I 4 

1=1 k=l v=l 
oo I 

E,AP) = EE a W M ££(p) , (3. 

1=1 k=l 



where N^Kp) and M^ l u (p) are obtained in 0(g 2 ) from the expansion of symmetrized 
k x / Wilson loops, 



w k > 1 = — - — y^w k ' l (x) 

VV W — AT3AT VV /i,J/W 

Jv cr Jv r x 



4iV 



-2M^(p)A«(p)A:(- P ) ) + Otf) . (3.9) 



Explicit expressions for N*r v (p) and M^'l(p) are given in the Appendix. 

With these definitions at hand, it is easy to evaluate the leading order expression 
for the energy density on finite lattices. Using the 0(g 2 ) expansion coefficients for 
k x I Wilson loops, given in the Appendix, one finds, 

= (tfW A-i + iV^A-i - 2M*J A"i) , (3. 10) 

and finally 

1 oo 2 



[i<u<4: 1=1 k=l 
3 oo I 



1 a CX) I 

6 ii=i i=i k=i 



The necessary inversion of the matrix A^ v and the evaluation of the sums ap- 
pearing in Eq. (|3.11| ) have been performed using Mathematica. We have explicitly 
checked that all our results are independent of the choice of the gauge fixing function 
<7 M (p) and the gauge fixing parameter £. In particular we have verified that in the 
infinite volume limit the 0(g 2 ) expansion coefficients are identical for all improved 
Wilson actions, 



lim lim S 1 '^ = lim lim S 1 '^ = - . (3.12) 

N T ^oo N a ^oo a A^7- — ► oo A^[j — > oo 2 

We also have evaluated separately the 0(g 2 ) expansion coefficients for the expecta- 
tion values of the different loops appearing in the improved actions. For S^ 1,2 ) our 



results on an infinite lattice agree with earlier calculations |TlJ, while we disagree for 
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S^ 2 ' 2 ) with results quoted in |]IEfl p[ For the expansion of the plaquette expectation 
value we find 



1 _ _ Re TrU^U^U^U+J 



N 



,N 2 -1 
AN 



+0(g 4 ) 



0.5 

0.366263 
0.392641 
0.358304 



Hl,l) 
I=(l,2) 
I=(2,2) 
I=(3,3) 
(3.13) 



3.2 Cut-off Effects in the Wilson Action 



2s ^, which, to 



For the one-plaquette action we use the gauge fixing function 
lowest order in a, is the momentum representation of the lattice version of d^A^. 
Choosing the Feynman gauge, £ = 1, we then obtain a diagonal propagator 



A^(p) = D(p)5^ with D{p)=AY J s 2 i 



(3.14) 



This leads to the familiar result for the energy density |10 



Si S A 



Dip) 



(3.15) 



Eq. (|3.15| ) is the starting point for a discussion of the cut-off dependence of 
the high temperature limit of the energy density. We will consider here the cut- 
off dependence in the thermodynamic limit (N a — ► oo) and will comment on the 
influence of a finite spatial volume later. In the thermodynamic limit the sum over 
the spatial momenta can be replaced by an integral. The remaining sum over p± can 
be performed explicitly [19]. One then finds, 



e 

rp4 



4(iV 2 - 1)N? 



(27T) 



2tt 



d 3 P 



VI + ^ 2 exp (N T x) 



(3.16) 



with x = 2 ln(u; + Vl + uj 2 ) and u 2 = sf + s\ + s\. The right hand side of the above 
equation can be expanded in powers of 1/N T = aT, which explicitly shows that the 

c The gluon propagator given in Ref.|l(| for the action S^ 2 ' 2 ) corresponds to a gauge, in which 
^y.,v = D^Sf!^. In general, this cannot be achieved with the covariant gauge fixing term given 
in Eq. (3.6) and therefore does require non-trivial ghost contributions already in leading order 
perturbation theory, i.e. for the ideal gas |18 . The ghost contribution has not been calculated in 
Rcf. [Ira. For this reason the 0(g 2 ) result quoted there is not the complete contribution to the 
expansion coefficient for the plaquette expectation value at that order in g 2 . 
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leading corrections to the continuum Stephan-Boltzmann law are indeed 0((aT) 2 ) 
T3 \. This expansion is somewhat tedious but straightforward. In next-to- leading 
order we fincQ 



T 4 T 4 J 15[ 63 N* 3iV 4 \N*)\ V ; 

In Figure |l]we show the cut-off dependence of e/T 4 obtained from Eq. (|3.16|) as well 
as the 0(N~ 2 ) and 0(N~ A ) results from Eq. (|3.17|) . As can bee seen the leading 




0.07 



1/N t 2 



Figure 1: High temperature limit of the energy density of an SU(N) gluon gas 
calculated with finite cut-off, aT = 1/N T , in units of the continuum limit result 
(solid curve). The dashed and dashed-dotted curves show this ratio calculated up 
to 0(N~ 2 ) and 0(iV~ 4 ), respectively. 

finite size correction term given in Eq. ( |3.17| ) accounts for more than 70% of the 
total correction already for N T > 6. 



3.3 Improved Actions 

Generalizing the discussion of the perturbative calculations with the one-plaquette 
Wilson action to the case of the improved actions is now straightforward, although 
the actual calculations turn out to be quite tedious because we did not succeed in 
finding a gauge fixing condition, which would lead to a simple diagonal propagator. 

d Note that the natural expansion parameter is irTa = n/N T , i.e. half the smallest, non- vanishing 
Matsubara mode. 
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Details of the perturbative calculations for the various improved actions considered 
by us are given in the Appendix. 



The generalization of the leading order perturbative expansion for the energy 



density, given in Eq. (|3.15) for the one-plaquette Wilson action, now reads 



T 4 2 V 



l)iV r 4 



R a (p)-R T (p) 
R d (p) 



(3.18) 



with functions R a , R T and R d defined in the Appendix for the different improved 
actions. 



We have evaluated the energy density from Eq. (|3.18|) numerically. Results for 
small values of N T obtained in the thermodynamic limit (N a — > oo) are summarized 
in Table |l|. We note that for all improved actions the deviation from the continuum 
results, ssb/T 4 ", is less than 1.2% already for N T = 6 which is reached with the 
Wilson action only for N T ~ 20. At this level of accuracy it is, in fact, important to 
consider also the influence of a finite spatial extent of the lattice. A calculation of 
e(N T , N a )/esB for finite N a indeed shows that the resulting infrared finite-size effects 
also are of the order of 1% for N a /N T ~ 4. Some results for finite spatial lattices 
are summarized in Table 0. 



Table 1: Deviations from the continuum ideal gas behaviour on spatially infinite 
lattices with temporal extent N T for the standard one-plaquette Wilson action and 
various improved actions. In the last row we give the coefficients of the leading 
cut-off correction term as defined in Eq. (3.19). The number in brackets gives the 
difference between the estimate for c 1 obtained for N T = 32 and the extrapolated 
value. 

e(N T )/e SB 



N T 


I=(l,l) 


I=(l,2) 


I=(2,2) 


I=(3,3) 


4 


1.495186 


0.986568 


1.087709 


1.044357 


6 


1.181566 


0.997528 


1.011994 


1.008095 


8 


1.086700 


1.000309 


1.003752 


1.001199 


10 


1.051708 


1.000253 


1.001582 


1.000320 


12 


1.034756 


1.000150 


1.000780 


1.000110 


c 1 


0.4762 


0.044 (2) 


0.178 (2) 


0.394 (19) 



The energy density calculated perturbatively from the improved actions will show 
a cut-off dependence starting at 0((aT) a ), i.e. corrections on lattices with temporal 
extent iV T will asymptotically be proportional to N~ a with a = 4 for the cases 



11 



Table 2: Deviations from the continuum ideal gas behaviour on finite spatial lattices 
with temporal extent iV T for the standard one-plaquette Wilson action and various 
improved actions. 





e(N T ,N a )/e SB 






N a = 


AN T 






N a = 


6N T 




N T 


HM) 


I=(l,2) 


I=(2,2) 


I=(3,3) 


I=(l,l) 


I=(l,2) 


I=(2,2) 


I=(3,3) 


4 


1.4833 


0.9747 


1.0758 


1.0325 


1.4917 


0.9830 


1.0842 


1.0408 


6 


1.1697 


0.9857 


1.0001 


0.9962 


1.1780 


0.9940 


1.0085 


1.0046 


8 


1.0748 


0.9884 


0.9919 


0.9893 


1.0832 


0.9968 


1.0002 


0.9977 


10 


1.0398 


0.9884 


0.9897 


0.9884 


1.0482 


0.9967 


0.9981 


0.9968 


12 


1.0229 


0.9883 


0.9889 


0.9882 


1.0312 


0.9966 


0.9973 


0.9966 



I = (1, 2) an I = (2, 2) and a = 6 for I = (3, 3) ; 



e(N T 



71 



2 r 



(N 2 - 1) 
T 4 v ; 15 



l + c 1 



7T 

K 



+ 0(N; 



(3.19) 



We have evaluated e/T 4 also for large values of N T , up to N T < 32, in order to 
verify this asymptotic behaviour. In Figure ||] we show estimates for the expansion 
coefficients, c 7 (iV r ), obtained at fixed values of N T by solving Eq. (|3.19| ) for c 1 . As 
can be seen, the cut-off dependence indeed shows the anticipated scaling behaviour. 
The coefficients c 1 have then been determined from an extrapolation to N T —>■ oo, 
taking into account the subleading correction term proportional to N~( a+2 \ These 
numbers are given in the last row of Table [[]. The quoted errors reflect the difference 
between the extrapolated value and the last calculated approximant for iV r = 32. 
We note that in the case of the improved actions S^ 1 ' 2 ^ and S^ 2 ' 2 -* not only the 0(N~ 2 ) 
corrections have been eliminated, but also the magnitude of the 0(N~ 4 ) coefficient 
has been reduced strongly when compared to the one-plaquette action. 

The improved actions lead to a drastic reduction of finite cut-off effects already 
for small values of iV r , where in the case of the one-plaquette Wilson actions the 
cut-off effects were not at all dominated by the leading 0(N~ 2 ) contribution, which 
has been eliminated in the improved actions. In fact, in the case of the improved 
action S^' 2 \ deviations from the continuum ideal gas behaviour are always less than 
1.5% already for iV T > 4, which is to be compared with the nearly 50% deviations 
at N T = 4 for the one-plaquette Wilson action. We note, however, that also in 
this case the leading N~ 4 dependence of the cut-off effects starts dominating only 
for iV T > 10. This is different for the actions S^ 2 ' 2 ^ and S^ 3 ' 3 \ where the leading 
correction dominates already for iV T > 6. 
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Figure 2: Estimate of the coefficient c 1 of the leading cut-off dependence, obtained 
by solving Eq. (3.19) for c 1 at fixed N T , with a = 2 for the one-plaquette Wilson 
action (triangle), a = 4 for / = (1,2) (crosses), / = (2,2) (circle) and a = 6 for 
I = (3,3) (squares). In the limit N T — > oo this quantity yields the coefficient cj 
given in Table 1. 

4 Monte Carlo Results for Improved Actions 




The perturbative calculations have shown that already the 0(a 2 ) improved actions 
I = (2, 2) and in particular I = (1, 2) provide a large reduction of finite cut-off effects 
in the high temperature (T — ► oo) limit even for small values of N T . We have used 
these actions as well as the 0(a 4 ) improved action, / = (3, 3), to analyze the cut-off 
dependence at non- vanishing values of g 2 , i.e. for T < oo. The recently performed 
simulations with the one-plaquette Wilson action on lattices with N T = 4, 6 and 8 



14fl do provide here a good basis for a comparison. 



In a first step we have evaluated the difference of action densities, studied per- 
turbatively in the previous sections, at finite values of the gauge coupling (3 on 
lattices of size 16 3 x 4 (24 3 x 4 in the case of the action S^ 2 ' 2 ^). We have performed 
calculations of the action differences for the SU(3) gauge theory, 

^) Q = 3(3N^(S I a )-(S I T }) , (4.1) 

at values of (3 = Q/g 2 well above the critical value for the deconfinement transition. 
Some results are given in Table |[ As can be seen, also at finite values of (5 the 
results for the different actions show a cut-off dependence, which is close to the result 
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calculated in the limit ft — > oo. In all cases the asymptotic value is approached from 
above. For the one-plaquette Wilson action this has been shown to be in accordance 
with the perturbatively calculated 0(g 2 ) correction to (e/T 4 ) ]nj. This behaviour 
should not be confused with the correction to the complete energy density, e/T 4 , 
which is negative. 

Table 3: Differences of space- and timelike action densities times iV 4 calculated 
on lattices of size 16 3 x 4 for / = (1,1), (1,2) and (3,3) and 24 3 x 4 in the case 
of / = (2,2). The last row gives the perturbative result in the limit ft — > oo 
(the corresponding value in the continuum is €sb/T 4 = 5.2638). Couplings have 
been chosen well above the critical coupling for the deconfmement transition, which 
occurs at ft c = 5.6908 (2) (I = (1,1), [20]), 4.0752 (13) (I = (1,2), [7]), 4.3995 (2) 
(I = (2,2), [this study]) and 3.7 < ft c < 4.2 (I = (3,3), [this study]). 



(e/T 4 ) 



ft 


(1,1) 


(1,2) 


(2,2) 


(3,3) 


6.0 


7.359 ( 12) 


5.706(57) 


6.513 (74) 


5.88 (12) 


10.0 


8.046 (102) 


5.42 (10) 


6.203 (90) 


5.79 (16) 


15.0 


7.837 ( 63) 


5.445(84) 


5.982 (79) 


5.58 (16) 


20.0 


7.806 ( 48) 


5.18 (10) 


5.860 (94) 


5.56 (13) 


oo 


7.808 


5.131 


5.707 


5.435 



A more detailed comparison with the one-plaquette Wilson action is possible 
through a complete calculation of the temperature dependence of bulk thermody- 
namic quantities along the line presented in |l4j . This also requires the calculation 



of the action densities on large, zero temperature lattices. As this is quite time con- 
suming we have performed such a complete analysis only for the action / = (2,2). 
On a lattice with four sites in the temporal direction this action leads, in the limit 
T — > oo, to finite cut-off effects which are of similar size as those found for the 
one-plaquette Wilson action on lattices with eight sites in the temporal direction. 
A direct comparison of results obtained with both action thus is possible. 

We have calculated the action densities (S^ 2 ' 2 ^)t and (S^ 2 ' 2 ^)q on lattices of size 
24 3 x 4 and 24 4 , respectively. An integration of their differences calculated at various 
values of ft yields the pressure^ 



y4 



" = K f"dft'((S^) -(S^) T ) , (4.2) 

0o J/3 



where fto is a value of the coupling constant below the phase transition point at 
which the pressure can savely be approximated by zero. In order to compare this 

c We refer to Refs. jH], for more details on the formalism. 
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calculation with corresponding results for the one-plaquette Wilson action, we have 
to define a common temperature scale. For this purpose we determine the relation 
between f3 and the cut-off via the 2-loop renormalization group equation, aA = 
R(Peg) using there an effective coupling constant obtained in terms of the action 
density 

N 2 -l 

a* = W% ■ (4 - 3) 

In the case of the one-plaquette Wilson action, this has been found to be a reasonable 



parametrization of the temperature scale [14]]. Combined with a determination of 



the critical coupling for the deconfinement transition on the 24 3 x 4 lattice, 

P e (N T = 4) = 4.3995 ± 0.0002 , (4.4) 

the temperature can be expressed in units of the critical temperature T/T c = 
R(P eS , c )/R(M- 

In Figure [3] we show the pressure calculated with the action S^ 2 ' 2 ^ on a 24 3 x 4 
lattice and compare this with the calculations performed with the one-plaquette 
Wilson action on lattices of size 32 3 x 4, 6 and 8 fT4|| . 

We note that at high temperatures, T>4T C , the pressure calculated from the 
improved action, S^ 2 ' 2 \ on a lattice with N T = 4 indeed is in good agreement with 
results obtained for the one-plaquette Wilson action on a lattice with N T = 8. This 
confirms that the improvement of the actions at g 2 = persists also at non-vanishing 
values of the gauge coupling. For temperatures closer to T c we find, however, results, 
which are compatible with those obtained with the one plaquette Wilson action. 
This is not too surprising. Already the analysis of the cut-off effects in calculations 
with the Wilson action has shown that the magnitude of the cut-off corrections does 
vary with temperature. 



5 Conclusions 



In the high temperature limit the thermodynamics of SU(N) gauge theories is dom- 
inated by high momentum modes. In lattice calculations these modes are most 
strongly affected by the finite ultraviolet cut-off. We have shown here that the use 
of tree level improved actions for the gauge fields indeed leads to a strong reduction 
of finite cut-off effects in thermodynamic observables like the energy density or the 
pressure. The perturbative analysis of the high temperature limit of the energy den- 
sity has shown that already on lattices with a temporal extent as small as N T = 4 
the cut-off effects can easily be reduced to a few percent. 
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Figure 3: The pressure in units of T 4 versus T/T c calculated with the improved 
action I = (2, 2) on a 24 3 x 4 lattice. This is compared with calculations using the 
standard one-plaquette Wilson action on lattices of size 32 3 x 4, 6 and 8 [14]. The 
dashed horizontal lines give the corresponding ideal gas limits for the (lxl) Wilson 
action and the solid line corresponds to the ideal gas limit for the (2, 2)-improved 
action (see Table 2). The arrow points at the continuum result for an ideal gas. 



A first Monte Carlo calculation of the pressure of an £77(3) gauge theory using 
an improved action has shown that this also removes the dominant cut-off effects 
at non-zero values of the gauge coupling. At high temperatures the calculations 
performed with an improved action on a lattice with N T = 4 lead to results which 
are compatible with those obtained with the standard one-plaquette Wilson action 
on a lattice with N T = 8, where this action leads to similarly small cut-off effects. 
However, closer to T c we also find that the inherently non-perturbative features of 
the deconfinement phase transition play a much more important role. Here the 
results obtained with a tree level improved action coincide with those obtained with 
the standard Wilson action on the same size lattice. Clearly in this temperature 
region there is need for further improvement of the lattice regularized Yang-Mills 
action. This can, for instance, be achieved by including one-loop corrections in the 
coefficients a^i of the improved action or by determining them through a Monte 
Carlo renormalization group approach at various values of g 2 . Currently we also 



investigate the properties of tadpole improved actions pi 



Of course, considerations similar to those made here for SU(N) gauge theories 
can also be applied to the fermion sector of QCD. In fact, in this case the influence of 
finite cut-off effects on the high temperature limit of thermodynamic observables are 
known to be even larger on lattices with small temporal extent. The calculation of 
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the equation of state of QCD thus should profit even more from the use of improved 
actions. An investigation of this claim is under way. 
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A Appendix 

We give here explicit expressions for the leading order weak coupling expansion of 
the action densities on finite lattices for the improved Wilson actions defined in 
Eqs. (^4] - p.4| ), respectively. For completeness we also repeat the corresponding 



results for the standard one-plaquette Wilson action. In all cases we have checked 
that the results are gauge invariant. 

In addition to the short hand notation s M = sin(p M /2) we introduce the abbrevi- 
ations 

M(p) = e*"-l , 

k— 1 

K k) (p) = • (a.i) 

1=0 

The expansion of a symmetrized k x I Wilson loop defined in Eq. (|2.2| ) to 0(g 2 ) can 
then be written as 

w k ' 1 = — - — rw y d) 

J V JV T X 

= 7T?E te(pMS(p)^(-P) + N^(p)A:(p)A a v (-p) 

p,a ^ 

rk.l /„\ a cl i 



with 



-2M* : l(p)A*(p)A a u (-p) + 0(g*) (A.2) 



K-» = \ (aW(p)aW(-p)a;(-Ip)a;(i p ) + {k~ i)) (A.3) 
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and 



Mft(p) = l( K e^-^ 2 A^\p)A;(lp)A^(-p)A^-kp) + (k <-> J)) • (A.4) 

We note that Mpl has non- vanishing matrix elements also for (// = v). In fact, 
due to the identity Ajf)(p)A~(Zp) = Af(p)A~(kp) one finds that M*'j, = iV*£. 
This term does not contribute in the expansion of Wilson loops. It is included in 
the definition of E^ v and cancels against the corresponding term in the diagonal part 
of the matrix G^ v = D^(p)5^ u — E^ u (p), which has been introduced in Eq. 
With the above relations we obtain for D^p) and E^ v (p), 



oo I i 

• //:. 



{=1 k=l v=\ 

n I , p fc-1 -| 4 

2 EE a M( fc + 2E(^-j)cos(j'PM) E sin2 (W 2 ) 



(=1 k=l 



3=1 



n i 

^(p) = EE fl *^(p) ( A - 5 ) 

i=l k=l 

n I ,k—l I— 1 

= 2 e E a Ki E sin (( 2 * + l )pJ 2 ) E sin ((2j + 1 W2) 

Z=l fc=l \=0 j=0 

+(Jfe <-> Z) 



Adding the gauge fixing functions ^ one can construct the propagator matrix 
A^j, (Eq. ( |3.6|) ), which we invert using Mathematica. 

The explicit results for the standard one-plaquette Wilson action have been dis- 
cussed in Section 3.1. We give here the result for the expansion coefficient, S 1 ^, 
for improved actions, which we write in the form 

01,(2) _ f R »AP) (A fi N 

b ^ ~ h R d ( P ) • { } 

In Eq. Q3.1S| ) we have introduced the functions R a and R T , which are defined as 

Ra- = g(-Ri,2 + -Ri,3 + -^2,3) 5 

Rr = 3(^1,4 + ^2,4 + ^3,4) • (A.7) 
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In order to present the results in a convenient form, we have chosen a specific 
form of the gauge fixing function 



I = (1, 1) : g„ = 2 ^ 

I = (1,2) : <fc = 2 (a„ + ? aj) 

I = (2,2) : g, = 2(s, + sl) 

1 16 

J = (3,3) : 9ll = 2(s lt + -s 3 fl + —s 5 fl ) 



(A.8) 



In the following we use the convention that all indices specified are mutually 
different. Permutations of indices that lead to identical expressions are counted only 
once. Whenever summations over permutations of indices are needed the relevant 
sums are specified explicitly. We present the result for the (//, v) component of the 
action. The remaining two directions of the 4-dimensional lattice are denoted by 
(i,j). In those cases where indices can take values from the set (fi, v) as well as (i,j) 
we denote them by k, I, m and n. With these conventions we find 



R 




4 



R d 



(A.9) 
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J = (1,2): 
R d 



3 (3 + sl + sl) 

DiD^sM + s v g v f + {D^l + D v sl)(D ig ) + D 3 g 2 ) 

Y,9 2 k DiD m D n 



(A.10) 



fc=i 



J=(2,2): 



R a 



I =(3,3): 



^ + 4 + 4-814) 

4 



k=l 



64 



(A.ll) 



128 



DiDjis^ + s u g u f + {Dpsl + D v s 2 v )(D ig 2 + D j9 2 ) 



256 



+ Y8^( D ^ + D » s » )Ai ' jBj > 

256 



18225 
4096 
'^225 



2 2/ 2 2\2f 2/ 2 2\/ 2 2\ , 2/ 2 2\ / 2 2\V 
- Si) ^(S M - SJ(S^ - Sj) + - Si ){s u - Sj ) j 



4 256 3 4 

J2 9lDiD m D n + E 



/»■ i 18225 fc=i /=fc+i 



40Q6 4 



225 



fc=i 



(A.12) 



Here we also have used the abbreviations 



Bk,i 



sl(sl-s*)(s 2 k s 2 + 3(s 2 k + s 2 )) 

4(4 - * 2 )(-405 + 120(4 + 4) + ^44) 



(A.13) 
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